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Abstract 

A (k,g)-graph is a fc-regular graph of girth g, and a (k, g)-cage is a (fc, g)-graph of mini- 
mum order. We show that a (3,ll)-graph of order 112 found by Balaban in 1973 is minimal 
and unique. We also show that the order of a (4,7)-cage is 67 and find one example. Fi- 
nally, we improve the lower bounds on the orders of (3,13)-cages and (3,14)-cages to 202 and 
260, respectively. The methods used were a combination of heuristic hill-climbing and an 
innovative backtrack search. 

AMS Subject Classifications: 05C25, 05C35 
Keywords: cage, regular graph, girth 

1 Introduction 

A k-regular graph is one in which every vertex has degree k. The girth of a graph is the length 
of a shortest cycle. A (fc, g)-graph is a regular graph of degree k and girth g, and a (k, g)-cage 
is a (k, g)-graph of minimum possible order. 

The problem of determining this order and identifying the corresponding cages has been 
extensively studied. See [3] for a current survey. The cases where the order of the cage was 
known precisely before this work can be summarized as follows. 
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1. Regular graphs of degree 3 for girths up to 12; 

2. Girth 5 graphs for degrees up to 7; 

3. Girth 6, 8 and 12 graphs for degree one more than a prime power; 

4. The case of degree 7 and girth 6. 

In this note, we add the case of degree 4 and girth 7 to the list, and also show that the 
(3,ll)-cage is unique. 

The order, but not the uniqueness, of the (3, ll)-cage was previously announced in [5]. 

2 Backtrack search 

The uniqueness of the (3,ll)-cage, and the lower bound of 67 for the (4,7)-cage, was proved using 
the program mentioned in [SJ but not described in detail there. We will provide that description 
here using the (4,7)-cage as an example. 

Consider the construction of 4-regular graphs of girth at least 7 and order n > 53. The 
vertices at distance at most 3 from some fixed vertex form a tree T with 53 vertices as in Figure [TJ 
(In the case of even girth we would root the tree at an edge rather than a vertex.) 




Figure 1: The initial tree T 

Let W be the set of n — 17 vertices consisting of the 36 leaves of T and the n — 53 vertices not 
in T. The task is thus to add additional edges within W so that the resulting graph is quartic 
and girth at least 7. We can do this in standard depth-first manner, starting with the tree and 
adding one edge at a time. For addition of a new edge vw to be valid, v and w must have degree 
less than 4 and be at distance at least 6; these properties are easily monitored. 

At each stage in the search, we choose one vertex of degree less than 4 and try all the 
possibilities for joining it to other vertices. We choose the vertex whose options are the most 
restricted, as experiments showed this heuristic to be a good one. 

Without additional improvements, this search is far too expensive due to multiple equivalent 
subcases. The most obvious source of equivalence is the automorphism group of the tree. Thus, 
the two choices a, b in Figure [5] are clearly equivalent and there is no need to try b as the first 
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choice once a has been tried. The simple structure of the tree allows examples of this type of 
equivalence to be monitored efficiently without explicit computation of automorphism groups. 



a 




Figure 2: Two equivalent choices 

However, when a moderate number of edges have been added, discovery of automorphisms is 
considerably more difficult. More importantly, there are equivalent subcases in the search that 
do not derive directly from automorphisms of any of the graphs that are constructed. Consider 
Figure [3l 




Figure 3: Complex equivalence 

Suppose we have completed the part of the search that begins with the three edges a, b, c in 
the upper diagram of the figure. This means that, up to isomorphism, we have already found all 
the quartic graphs that contain the upper diagram as a subgraph. Now suppose that at a later 
point in the search we have added the edges a, b, d, e, f shown in the lower diagram of Figure [3j 
and are about to consider further edges. If we add edge g, then we will not find any new quartic 
graphs, since T + {e, /, g} is isomorphic to T + {a, b, c}. Therefore, we can mark g as ineligible. 
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This applies whether we consider g as the very next edge to add or we consider adding it to 
some non-trivial extension of a, b, d, e, f. For efficiency reasons we separate these two cases. 

We now define the pruning process formally. Each node of the search tree corresponds to a 
graph T + E, where E is a set of edges within W. The subsearch rooted at T + E is the subtree 
of the search tree rooted at T + E, and T + E is completed when we have finished the scan of 
that subtree. 

Pruning Rule. Suppose that T + E is a search node and vw is an available edge (that is, v 
and w have degree less than 4 and are at distance at least 6 from each other). Suppose 
further that there is a subset E' C E + vw such that T + E' is isomorphic to a node which 
is already completed. Then the subsearch rooted at T + E can avoid adding the edge vw. 

The pruning rule is applied using the nauty graph isomorphism software [3]. Since subgraph 
isomorphism testing is required, the rule is very expensive compared to the time otherwise 
required to process one node of the search. In practice we limit it to only certain E and certain 
v,w. Specifically, we define two levels t\ > £2 > 0. For search nodes T + E such that \E\ < £2, 
we apply the rule for all v,w. For li < \E\ < £\, we apply it only to those v,w such that vw is 
a candidate for the very next edge to add. 

For larger £i,£%, the number of nodes in the search tree is reduced but the time expended in 
pruning the tree is increased. So there is some optimal compromise. 

A further technique is needed in very difficult cases (such as n = 66) , for which it is desirable 
to divide the computation across multiple processes (perhaps on different computers). The 
division is accomplished in the usual fashion: for some level £q, treat the subsearchs rooted at 
T + E for \E\ = £q as independent. Subsearch i is assigned to process j if i = j (mod N), where 
N is the number of processes. 

Each process involved in the search computes the whole search tree to level £q and only its 
own subsearches at higher levels. Since the application of the pruning rule is hard to achieve 
across multiple processes, we use £q > £\. Higher £q also tends to share the load more evenly 
between processes. These considerations, however, imply that a large part of the tree, including 
all the rule applications, are repeated in every process, which would seem to prevent use of 
large £1,^2 due to the expensive nature of the pruning rule. The solution is to conduct the 
computation in two phases. In the first phase, executed on a single processor, the search tree 
is computed up to level l\ while applying the pruning rule. During this phase, the results of 
all the rule applications are recorded in an audit file. Then, in the second phase where the full 
search tree is computed, the pruning rule is applied at almost no cost by following the audit file. 
If Iq is not too high, the cost of the division of the second phase into parts is negligible. 

For n = 66, we used £q = 22, i\ = 16 and £2 = 6. The total number of nodes in the search 
tree was 318,904,129,273,923. Using a large mix of computers ranging from 650 MHz to 3 GHz, 
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the total time was 96 years (105,000 nodes per second). This cost was about 5 times the cost 
for n = 65. 

Improvements to the program after the computation finished resulted in a 4- fold speedup. 
Nevertheless, the computation of all the cages of order 67 will not be feasible in the near future. 

Exhaustive computation of the (3, ll)-cages, showing that the graph found by Balaban [lj 
is the only one, required 17 years on computers averaging about 300 MHz. 

With lesser expenditure of cpu time, we also showed that a (3, 13)-cage has at least 202 
vertices (improved from 196), and that a (3, 14)-cage has at least 260 vertices (improved from 
256). The first result, and the bound 258 for the second, were previously announced in [5]. 

To emphasise the efficiency of this method, we compared it against the previously best code 
for cubic graphs [2]. For girth 9 and 58 vertices, our approach was 52 times faster, while for 
girth 11 and 104 vertices our approach was 856 times faster. All of the results we highlight in 
this paper were previously infeasible. 

3 A Graph of Degree 4 and Girth 7 

As noted in the previous section, we showed exhaustively that there are no (4, 7)-graphs with 
66 or fewer vertices. To complete the proof that (4, 7)-cages have 67 vertices, we present an 
example. 

Our example on 67 vertices is given below as an adjacency list. The graph has an automor- 
phism group of order four, isomorphic to Z2 x Z^. In the action of this group on the vertices, 
there are 11 orbits of length 4, 11 orbits of length 2, and one fixed point. 

The graph was constructed using a hill-climbing algorithm that begins with an empty graph 
on 67 vertices and adds edges, one at a time, while not violating the degree and girth conditions. 
An outline of the algorithm follows. 

while there are vertices with degree < 4: 

while there are edges that can be added: 
for each edge that can be added 

compute the degree sum of its vertices 
pick the edge with the largest sum 
choose randomly in case of ties 
add the winning edge 
if all vertices have degree 4: 
save graph 
exit 

delete 1, 2 or 3 random edges 
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Figure 4: Adjacency list for a (4,7)-graph of order 67. 



A few of the steps require some explanation. 

1. Edges are added between vertices with maximum possible degree sum. Doing this in the 
early stages of the algorithm appears to be essential to finding a cage. We ran our program 
many times with this condition removed, without success. 

2. When edges are added, a record is made of when they were added. Time is measured in 
trips through the outer loop. 

3. When edges are chosen for deletion, it is done in one of two ways. Either the probability 
that an edge is chosen is proportional to its age, or it is inversely proportional to its age. 
These two modes are alternated, each being used for a few thousand time periods. Very 
recently deleted edges are not chosen for reinsertion. 

A program implementing this method has been run several hundred times. It succeeds in 
finding a graph approximately twenty percent of the time. Each time it has succeeded, it has 
found the same graph. 



6 



References 



[1] A. T. Balaban, Trivalent graphs of girth nine and eleven and relationships among cages, Rev. 
Roumaine Math., 18 (1973) 1033-1043. 

[2] G. Brinkmann, B. D. McKay and C. Saager, The smallest cubic graphs of girth nine, Combin. 
Probab. Cornput., 5 (1995) 1-13. 

[3] G. Exoo and R. Jajcay, Dynamic cage survey, Electron. J. Combin., 2008, 48pp. Available 
at 



http : //www. combinatorics . org/Surveys/dsl6 .pdf 



[4] B. D. McKay, Practical graph isomorphism, 10th. Manitoba Conference on Numerical Math- 
ematics and Computing (Winnipeg, 1980), Congr. Numer., 30 (1981) 45-87. Software at 
http : //cs . anu . edu . au/ ~bdm/ nauty/ . 



[5] B. D. McKay, W. Myrvold and J. Nadon, Fast backtracking principles applied to find new 
cages, 9th Annual ACM- SI AM Symposium on Discrete Algorithms, January 1998, 188-191. 



7 



